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Abstract 

It was shown recently that Frobenius reduction of the matrix fields reveals interesting 
relations among the nonlinear Partial Differential Equations (PDEs) integrable by the 
Inverse Spectral Transform Method (S-integrable PDEs), linearizable by the Hopf-Cole 
substitution (C-integrable PDEs) and integrable by the method of characteristics (Ch- 
integrable PDEs). However, only two classes of S'-integrable PDEs have been involved: 
soliton equations like Korteweg-de Vries, Nonlinear Shrodinger, Kadomtsev-Petviashvili 
and Davey-Stewartson equations, and GL(N, C) Self-dual type PDEs, like Yang-Mills 
equation. In this paper we consider the simple five-dimensional nonlinear PDE from an- 
other class of S'-integrable PDEs, namely, scalar nonlinear PDE which is commutativity 
condition of the pair of vector fields. We show its origin from the (l+l)-dimensional 
hierarchy of Ch- integrable PDEs after certain composition of Frobenius type and differ- 
ential reductions imposed on the matrix fields. Matrix generalization of the above scalar 
nonlinear PDE will be derived as well. 



1 Introduction 

We continue study of the relationship among nonlinear integrable Partial Differential Equa- 
tions (PDEs) with different integrability algorithms in the spirit of refs. [H [2] , where the block 
Frobenious Reduction (FR) of the matrix fields has been introduced for this purpose. First of 
all, we recall that there are several classes of completely integrable nonlinear PDEs depending 
on their integration algorithm. Most famous classes are following. 

1. Nonlinear PDEs linearizable by some direct substitution (or C-integrable PDEs) [3[ HI 
El El d IE] - Mostly remarkable nonlinear PDEs of this type are PDEs linearizable by the 
Hopf-Cole substitution [9 J and by its multidimensional generalization [TO] . 

2. Scalar equations integrable by the method of characteristics [H] (or Ch- integrable PDEs) 
and their matrix generalizations [T2l [T3] . 

3. Nonlinear PDEs integrable by the inverse spectral transform method [HI [151 EH] an d 
by the dressing method [HI [TBI EEl EQ[ [21] (S'-integrable PDEs). We underline three 
subclasses of these equations: 

(a) Soliton equations in (l+l)-dimensions, such as Korteweg-de Vries (KdV) [HI [22] and 
Nonlinear Shrodinger (NLS) [23] equations , and soliton (2+l)-dimensional equa- 
tions, such as Kadomtsev-Petwiashvili (KP) [24J and Davey-Stewartson (DS) [25] 
equations. Hereafter we call this subclass SVintegrable PDEs. 
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(b) Self-dual type PDEs having instanton solutions, like Self-dual Yang-Mills equation 
(SDYM) [26l [27] and its multidimensional generalizations. This subclass will be 
called S^-integrable PDEs. 

(c) PDEs assotiated with commutativity of vector fields [23 1231 EDI ED E21 E31 EI!- 
These equations may be either scalar or vector (S^-integrable PDEs). 

It was shown in [2] that C and C/i-integrable matrix PDEs supplemented by the Frobenius 
reduction of the matrix fields lead to the proper class of S'-integrable PDEs. Thus, matrix PDEs 
linearizable by the Hopf-Cole substitution generate (2+l)-dimensional Si-integrable PDEs; 
matrix equations integrable by the method of characteristics generate ^-integrable PDEs or, 
in the particular case of constant characteristics, (l+l)-dimensional Si-integrable PDEs. 

However, the technique developed in [2] does not allow one to involve ^-integrable PDEs 
into consideration. In this paper we consider the particular example of Ss-integrable PDEs 
which may be incorporated into the above mentioned algorithm with minor modifications. 
Namely, we consider the following representative of S'a-integrable PDEs [34|: 

l^ziti ^Ziti r U j Z2^z\x U Zl U z ^ x 0, (1) 

whose Lax pair reads: 

V^(A;f) + \ijj Zi (\;x) + u Zi (x)ip x (X; x) = 0, 2 = 1,2, (2) 

where if) is a scalar spectral function, A is a spectral parameter and x = (x, Zi,Z2,ti, t 2 ) is a set 
of all independent variables of nonlinear PDE. We show (for the restricted manifold of solutions 
to the eq.flU)) that eq.flU) originates from the S^-integrable GL(N S ' 2 ,C) SDYM (iV 52 is some 
positive integer) with independent variables (zi, Z2, h, £2) after two reductions imposed on the 
matrix field: 

1. the Differential Reduction introducing one more variable x (DR(x)) 

2. the Frobenius Type Reduction (FTR) (or, in particular, just Frobenius reduction). 

Matrix generalization of eq.([T]) will be introduced as well. Remember that ^-integrable 
GL(N S ' 2 , C) SDYM originates from the (1+1) -dimensional C/i-integrable hierarchy of nonlinear 
PDEs through the Frobenius reduction [2]. All in all, the diagram in Fig.(l) illustrates the chain 
of transformations leading from the (l+l)-dimensional matrix N x N (where = 2NonoM) 
Ch- integrable hierarchy to the scalar S^-integrable PDE (1). 
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Fig.l The chain of transformations from the (l+l)-dimensional C/i-integrable hierarchy to the five- 
dimensional S 3 integrable PDE (1). Here N = 2N n M, N S2 = 2n M, I N s 2 and I M are N S2 - and 
M-dimensional identity matrices respectively, W is N x N matrix, are 2uqM x 2uqM matrices 
and are UqM x uqM matrices, v^O and are M x M matrices, and are defined 
by eqs.fllB 



Note that the reduction 

zi = x, t 2 = t\ = -Z2 = y (3) 
reduces eq.flT]) into the following one 

U x t ~t~ ^yj/ ~t~ ^x^xy 'Uy'U'xx 0> (4) 

which has been studied in [351 E3 E3 [38] . 

The structure of this paper is following. In Secj2]we recall one of the results of [2]. Namely, 
we describe a version of the dressing method relating the (l+l)-dimensional Ch- integrable 
hierarchy of nonlinear PDEs with the four-dimensional ^-integrable GL(N 32 , C) SDYM. The 
later admits differential reduction of the special type introduced in Secj3j After that, using FTR 
we derive assotiated five- dimensional matrix system of two matrix nonlinear PDEs in SecJH 
Scalar case of this equation yelds the ^-integrable eq.fll]). Solution spaces to the nonlinear 
PDEs derived in SecsJH and @] will be considered in SecEJ Conclusions are given in SecEJ 

2 Derivation of GL(N S2 ,C) SDYM from the hierarchy of 
matrix Ch- integrable nonlinear PDE. 

In this section we describe briefly the algorithm relating the simplest (l+l)-dimensional hierar- 
chy of the matrix N x N Ch- integrable nonlinear PDEs from one side and GL(N S2 , C) SDYM 
from another side. Here N = NqN S2 and Nq is an arbitrary positive integer. We start with the 
following linear equation [H [2] 

X A = W X (5) 
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with x given as a solution to the following linear PDEs: 

Xt n + Xz n A = 0, n = l,2,.... 



(6) 



Here x an d W are 2NqtiqM x 2NqHqM matrix functions and A is a diagonal 2NqUqM x 2NqtiqM 
matrix function. Parameters M, no and Nq are arbitrary positive integers. Meaning of these 
parameters is clarified in Fig.l. Eq.([SJ) must be viewed as algebraic equation for W with given 
X and A: 



W = xAx' 1 - 



(7) 



Note that the diagonal matrix A is not constant in general. In fact, compatibility condition of 
the eqs.((5]) and © reads 



x (A tn +A Zn A) = (W tn + W Zn W)x, n 
which suggests us the following nonlinear PDEs for A 

At, 

and for W 

W tn + W Zn W = 0, n=l,2,.... 



1,2, 



A 2n A = 0, n 



1,2, 



(8) 



(9) 



(10) 



Eq. © is integrable by the method of characteristics [121 E 12] • Now both x an d A are fixed, 
so that W (solution to eq. (|T0l ) may be found using eq.©. 

Similar to [2], in order to derive SDYM we take W in the Frobenius form: 



W 



( 

l2Mn ^2Mn 

§2Mno hMno 

\ §2Mn §2Mno 



W (No-l) W (N ) ^ 
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;n) 



where Ij and 0j are J x J identity and zero matrices respectively, are 2uqM x 2n M 
matrix functions. Let matrix parameter A be given in the following block-diagonal form 



A = diag(A (1) 



A (N 0)] 



(12) 



where A^ are 2hqM x 2uqM diagonal matrices. Substituting f TTTT) and ( 1T21) into eq. (JHJ) we 
obtain that \ must have the following block structure: 



/ 



X 



X 



(i) 



x (i) (A (i))-i 
V x (i)(A( 1 ))- Ar (»+ 1 



X 



(No) 



X 



(No) 



( A (iVo 



X 



(No) 



( A (iV 



)\~N +1 



(13) 



where ar e 2n M x 2n M matrix functions. In turn, 2NqUqM x 2NqUqM matrix equation 
(151) reduces to the following set of 2rioM x 2rioM matrix equations: 



x o-) A (i) = ^w«x w (A £, ' ) )- <+1 , J = l,...,iVo, 



(14) 
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while eqs.fjS]) and (J9]) yield 

x% + xS aU) = 0, ( 15 ) 

Ajf +Ag)A^ = 0, j = l,...,A , n = l,2,.... (16) 

Then compatibility condition of eqs. (114j) and (f!5l) yields the following chain [U [2] : 

W S + Wj£>W® + W {i+l) = 0, z = l,...,A , ^ (JVo+1) = 0, n = l,2,..., (17) 

which may be obtained directly after substitution eq. liTTl) into eq. fTTOj) . Putting i = 1 and 
eliminating using two equations (ITT)) with n = 1,2 one derives GL(N S2 , C) SDYM, A 52 = 
2n M: 

Wfl - W?\ + W 2 (1) W 2 (1) - ^j 1 ^ 1 ) = 0. (18) 

3 Differential reduction of \ and assotiated system of 
nonlinear PDEs 

Introduce one more reduction. Namely, let matrices x an d ^ have the following block 
structures: 

X U) = { jw-i) ^( 2i) J , A 0) = diag(A (2i_1) , A (2j) ), j = 1, . . . , N , (19) 

where \p( m ) are noMxnoM matrix functions and AM are uqMxuqM diagonal matrix functions. 
Eqs. (USD and (JTSD yield: 

+ ^ ] A {m) = 0, (20) 

A£° + A^A^ = 0, m = l,...,2A , n = l,2,.... (21) 

x-dependence of \I/M is introduced by the following second order PDE: 

£(°) : = $M = a ^ (m) A (m) + i^ m) + //^ (m) , m = 1, . . . , 2A , (22) 

where a, z/ and // are noM x noM diagonal matrix parameters. These parameters must be 
constant and A*™) must be independent on x. In fact, compatibility condition of (120)) and (122)) 
yields, along with (I2~T1) . the following conditions: 

AM =0, m = 1,...,2A , a, i/, /i are constant diagonal matrices, m — 1, . . . , A . (23) 

The block structures of x^ an d A^') ffl9|) suggest us the relevant block structures of W^': 

^ <24) 

where q^' and pV' are noM x noM matrix functions. Now set of 2uqM x 2uqM 

eqs.t lll)) may be written as two n M x 2N n M equations: 

iVo 

: = = ^ (v (i) ^ x + w (i) ^A- i+1 , (25) 
i=i 

N 

£ (2) := ^ X A = (q®V x + A" i+1 , (26) 

i=l 
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where 

* = -,* (2Aro) )- (27) 

Compatibility condition of eqs. fl25l) and fl26|) . 

^ 1} =^ (2) , (28) 

yields the expressions for p^> and g«) in terms of and w^'\ 

p(?) = + w Cj)j/ + ^ +1 ) a + dWoicW, g w = + + w U) + v {1) av {j) . (29) 

Thus only two blocks of W^' are independent, i.e. w^' and v^'. 

Matrix equation (125]) may be considered as the uniquely solvable system of scalar linear 
algebraic equations for the functions v® and w®, i = 1, . . . , N , while eq.(l2"oT) is the consequence 
of eq. fl25p . In order to derive the nonlinear PDEs for v® and w®, we turn to the compatibility 
condition of the eqs. fl20l and fl25l) : 

No N 

S^+S^A =► ^ ^)^A- i+1 + ^ E( 0i ^A- i+1 = 0, n = l,2, 

i=l i=l 

which generates the following chains of nonlinear PDEs for v® and w®, i — 1, . . . 

£f ) : = + + + A + uWavW) + w «„(0 + ug+i) = o, 

v (No+l) = Q 

M«) : = w « + ^W*) + + v awo + t,(*+l) a ) + + - 

lb L n z n \ jj • / An 6-n 

^o+i) = 0j n = i ;2 ,.... 



(30) 

,Wo: 

(31) 

= 0, (32) 



In addition, we must take into account the compatibility condition of the eqs. fT2"2"|) and (T2"5]) 

No 

££> = £ (0) A ^{E^ x + ^( w )*)A- i+1 = 0, (33) 

8=0 

which gives us the non-evolutionary part of the system of nonlinear chains, % = 1, . . . , Nq: 

EM := v® + [w®,u] + 2w® + 2vfv - wf + [v {i) , fj] + [v®,u)v + [v {1 \ u]av {i) + (34) 
2vWav® + [w^\ a]v® + [v {1 \ a](w® + v® + v®u + v {l) av®) + [v (i+1 \ a] = 0, 

E^ := w® + [w®,fj] + [w w ,a}w (i) +2v®/i+ [v®,v]fi + [v {1) , u}aw {i) + (35) 
2v^aw {i) - vwf + [v {1) , a](w® + v® fi + v^aw® + v {i+1) a) + a] + 

2^ +1) a+ [v {i+1 \v]a = 

Of course, the system (I31f32|34|[35|) may be derived directly from the system (fT7|) using reduction 
(J24]). However, the derivation of this system as the compatibility condition of the linear system 
(I20|I22|I25I) is more illustrative. 

The complete system of nonlinear PDEs for v® and w®, i = 1,2, is represented by the 
eqs.f l31|l32l with fixed n (say n — 1) and i — 1 and by the eqs. f)34|35p with i = 1: 

Ei n) := + «J)(u)W + 4 1 ) + v^u + vMavM) + w^v^ + v&> = 0, (36) 
E [ 01) := wg> + v^(w^ + v^/jl + V V>aw<U + v<®a) + h,J# + w&> = 0, (37) 
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EK 11 ) : = V W + [ w W,u]+2w£ ) +2vMv-uvj» + (38) 
[v {1 \v](v + av {l) ) + 2v ( £ ) av {l) + 

a]w (1) + a](w (1) + + v®v + v {1) av {1) ) + [v {2 \a] = 0, 

f( 01 ) : = u, w + ; M ] + [ w w } a }w w + 2^V + , Ha* + > + ( 39 ) 

2 V ( 1 W 1) - ^« + a](«;W + v (1 V + f (1, aw (1) + v {2) a) + [w/ 2) , a] + 
2t/ 2) a + [v (2) ,iy]a = 

The scalar version (uqM = 1) of the system (I36ti39[) reads: 

E 1 = +vW(wM +vW + wW +av(%W)+w (40) 

+ "gW + + av^wW + av<®) + tugV 1 ) + w$> = 0, (41) 

u£ + 2w« + + 2av^v {1) = 0, (42) 

u>S2 + 2 ^V + 2av^w {1) - vw® + 2av {2) = (43) 



4 Frobenius type reduction and assotiated higher dimen- 
sional systems of nonlinear PDEs 

It is remarkable that the chains of nonlinear PDEs ( l3~H32ll34f35l) admit the following FTR: 

v® = {v ( ?> kt >, k,l = l,...,n }, w {i) = {w {i ' H \ k,l = l,...,n }, (44) 
V (W) = 6 klV M) + 5 k{ i +ni{i)) I M , 

W {i ' H) = S kl W {i '' l) + 4(/+n 2 (i))^M, 

where n is a positive integer parameter, n\(i) and n 2 (i) are arbitrary positive integer functions 
of positive integer argument, and ti)" ! ' are M x M matrix fields. In particular, if rii(i) = 
ri2{%) = 1, then this reduction becomes Frobenius one |2], which is shown in Fig.l. Eq. (l44|) 
requires the following diagonal block-structures for a, v and [i: 



v = di&g(p,...,i>), fi = diag(//, . . . ,fi), (45) 





(u, \l and a are M x M diagonal constant matrices) and the appropriate block-structures for 
the functions and A™: 

^(0 = {$(H i n , m = 1, . . . , n }, (46) 
A w = diag(A (/;1) ,...,A (Z;no) ), / = 1,...,2JV . 

Here \I/(' ;nm ) are M x M matrix functions and A^ ;n ^ are M x M diagonal matrix functions. 
Eq. ([25]) reduces to the following one: 

* (,;nro) A^ = [( S niv {i ^ + S n(j+nm )^ m) + (47) 

i=i j=i 

/ = !,.. . ,2A^ n,m = 1, . . . ,n Q , 
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while eqs. ( |20|22l) and ( 12T1I231) yield respectively 

*ff' = a¥ l ' jm) A {m) + v¥^ m) + /i^ (/;jm) , (49) 
A£ m) + A^ m) A (/;m) = o, A^ m) = 0, (50) 
n=l,2, l = l,...,2N , j,m = 1, . . . ,n . 

Then the chains of nonlinear PDEs (I31|32II34|[35|) get the following block structures: 

Eg* = {EW»6 kl , k,l=l,...,n } = 0, (51) 
E(mi) = { E(mi;i) 5ku k ,l=l,...,n } = 0, 

m — 0, 1, z = 1, . . . , AT . 

where 

:= + v m( w M + + yfeQj?) + (^V 1;1 ) + ^PDjS^M + (52) 
4f^ (1;i+niW) S + ^%^ + (Qf ;0 ),„ = 0, 
£;(oi;0 := w g0 + v (i;i)( w (^) + ^ + v m;Da) + ( v g%M) + vif +n ^)aw^ + (53) 

E (ii;l) := „(*;/) + p] + 2«;(< i0 + 2^P - i/i;<W + /i] + [v {i;l \ u]u+ (54) 

[v^.j/Jov^ + [v {1 ' l+ni{l) \v]a + 2v£' 1) av {i ' l) + 2v£' l+ni(i)) a + [w^, a)v^ + 
[vW\a](wM + vM+vMu) + ([v^\a]v^ + [v^ 1+n ^\ ~a])~av^ + 
[v^\a]v^ +ni ^a+[Q^ t) ,a] = 0, 

E<f*$ := W W + [u;< i; V] + (55) 
^(i ; j+n 2 (i)) } ~]~ + 2v^ l) aw^ l) + 2v£' l+n2(i)) a - vwf l) + 

[v^ 1 \a]v^ l+n ^h + 2v^~a+ [v^ l \v\a+ [Q% l \~a] = 0, 

g(<;0 = ^(y+niB+mfl))- _j_ w (l;I+»»a(<)) + v (l;'+ni(i))p + ^CM+niW) + v (*+i;Q (56) 
gW) _ y(l;J+na(i)+m(l))5 + y(l;J+ni(i))^ _|_ ^CM+ni^+l))^ + ^(ljH-Mi)) _|_ w (*'+ 1 :0, 



We see that after the reduction ff44l) the chains of nonlinear PDEs (I31ll32|34|[35l) acquire one 
more discrete variable. 

To write the complete system of nonlinear PDEs we, first of all, put i — I — 1 in the chains 
of PDEs (I52H55P and take two values of n (n — 1,2) in eqs. ()52|l53l) (remember, that deriving 
the complete system (I36H39I) in Secj3]we fixed n = 1 in eqs. (13111321) ) : 

Si 11 ' 1 ) := vff + 4fV 1;1) + + ^u) + (v^vW + (57) 
u (i ; i+n a (i)))^(l;D + ^ryia+mCi))- + ^^(i;!) + (Q^ 1 ))^ = , 

EgW := iu£ ;1) + <4i ;1) («4 1;1) + + ^ (2;1) a) + (t£ ! V l!l) + (58) 

v^ 1+n ^)aw^ + v^v^+^a + w^w™ + (Q? ;1) k = 0, 
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MP 1 ** := + [w^\D] + 2wW + 2v^V - Vv^ + [v^\p\ + [v^,u]u + (59) 



v 



(i;D p] 5i; (i;i) + ^fti+md))^^ + a^W 1 *) + 2^ 1+ ™ 1 ( 1 »a + 



([v^a]^ + [^ 1+ni «),a])a^ 1;1 ) + [^ 1;1 ),a]^ 1;1+ni «)a+ [gS 1;1) ,S] = 0, 
E<f>W : = w (i/) + [ w a ; i)^] + [ w (i;i),a]^ 1 ) + 2^ 1 )/i+ (60) 

+ [u^ 1 ), i>]aw (1;1 ) + [ v (i;i+« 9 (i)) j + 

([r^U]^ 1 ) + [^( 1;1+ni ^,a])a«;( 1;1 ) + a]v^ 1+n ^a + 2vf^a + 

[v<* 1 \P]a+[Q§* 1) ,a] = 0, 



where 

Q (W = u (lil+an 1 (l)) 5 + w (l;l+n 2 (l)) + v (l;l+n l( l)) ~ + w (l;l+n l( l)) + u (2;l) j 
g(l;D = v (l.,l +ni (l) + n 2 (l))~ + v (l;l + n l( l))~ + w (l;l+n l( 2)) ~ + w (l;l + n 2 (l)) + ^(2;!)^ 

One can eliminate Qj 1 ' 1 ^ from the system (I57II59I) resulting in the following complete system of 
two PDEs for the matrix fields u = and q = w^ 1;1 ' + v^ l ' 1+ni ^ 1 '' ; 

(E[ 1V >\ 2 - {E^\, = 0, (62) 
[E[ lv ' 1 \~a}-(E^) Zl =0. (63) 

Three equations (I58ll60p are not important because they introduce three more fields (w^ 1 ' 1 ', 
Qz^ and p = i/ 2;1 ) + ^C^+M 1 ))) which do not appear in the system fl62|63l) . In particular, if 
a is a scalar, eq. (l63l) must be replaced by the following one: 

= 0. (64) 



Consider the scalar case (M = 1). Then only field u = v^ 1 ' 1 ' remains in the eq. 
which becomes eq.fll]). Eq. (163]) is not important in this case. Thus the system ( 16211631 may be 
considered as a matrix generalization of eq.([T]). Its integrability must be studied more carefully. 

All in all, we have the following chain of transformations relating four systems of nonlinear 
PDEs having different integrability properties (compare with Fig.l): 



eq . m cq JP eq . $® *M eqs. (JSH1 -ffl eqs. (mm scal ^ ase eq. © . (65) 



5 Solutions to the nonlinear PDEs 

5.1 Solutions to the system of nonlinear PDEs (I36H39I) 

Solutions to the system of matrix nonlinear PDEs (I36ti39[) may be written in terms of the 
functions \I>( m ) and M m \ m = 1, . . . , 2N , taken as solutions to the linear systems (I20II22P and 
(121f23p respectively with n — 1. They read: 

E^^i-A^Oe^^, (66) 

i=l 

Ef{ Zl -kfti), a,P = l,...,n M, I = 1, . . . , 2N . (67) 



A 
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where tp^iy) and Ep\y) are arbitrary scalar functions of single scalar variable, x = (x,z\,ti) 

is the list of all independent variables of the nonlinear PDEs, k^g are the roots of the charac- 
teristics equation assotiated with eq. fl22l : 



k 2 - a a kf - v a k - n a = (68) 



*$» = 2 + V (^) 2 + 4flA /3 + 4 ^ J > k ah = 2 - V ^ + 4aA /^ + 4 ^ 

Now, we can use eq. (1251) in order to find and u^ 4 -*. The simplest nontrivial case corre- 
sponds to N = 2. Then eq. (j2"5"]) reduces to the following four matrix equations: 

4 —i+l 
^M^(m) = ^W^M + w (i)^(m)^j ^(m) j 2 ; m = ^ 2 , 3, 4. (69) 

8=1 

These equations, in general, are uniquely solvable for the matrix fields and w;®, % = 1,2. 
We do not represent the explicite expressions for them. 

5.1.1 Dimensionality of the solution space 

To define restrictions, generated by our algorithm, on the solution space to eqs. (I36H39P we, first 
of all, write equation (1251 in the following compact form: 

yWjtfm) = ^yWn(H ) m=l,...,2JV , (70) 

71=1 

where = (v^ w^) and IT mnj is the following 2N n M x 2N n M invertible operator: 

(-pr(nm) \ 
^ nm) J , rt nm *> = *Wa(-^i), n = 1, . . . , N , m = l,...,2N . (71) 

Introduce operator IT" 1 -?) by the formula 

2N ° 

n(-)n^') = 6 nj i noM , n<™> = (n< mn) nf n) ), n , j = 1, . . . , tv . (72) 

m=l 

Then eq. ljTO]) yields 

v (n) = $HAHfjM ) w (n) = ^NAHfl^, n = 1, . . . , 7V . (73) 

m=l m=l 

Now let us turn to the eq. (1661) . We represent arbitrary functions ip^g{y) in the following 
integral form: 

.(M)/„,\ _ / j„j,(W 



€^(?/)= dq^>;>(q)e™, i = l,2, / = 1, . . . , 2iV , a, /3 = 1, . . . , n M (74) 
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where t/t^ , (g) are arbitrary functions of argument and one integrates over the whole space of 
the parameter q which is complex in general. Then eq. (!66|) may be written in the following 
form: 

2 - 

= E d ^«£(vy Zl ~~ A ^ 1+k ^ X i i = l,...,2JVb, «^l,..,n M. (75) 

i=l ^ 

Substituting expression (17511 for into eqs. (173l) we write them as follows: 

n M 2N 2 „ 

= E E E / ^^(gJAHe^-M"^^-^^)^, (76) 

7=1 m=l 8=1 
n M 27V 2 „ 

= EEE / w^(?)A^ e ^-^ m) *+^°«(n^ 

7=1 m=l i=l 

j' = l,...,JV , a,/3 = 1, . . . ,n M. 

These formulae might be considered as integral representations of and if fin (n = 1,2) 
would not depend on x. However, IT" 1 -^ do depend on x, so that eq. (IT6I) has more complicated 
sense. Nevertheless, in order to estimate the dimensionality of the solution space we consider 
eqs.f )76p as integral representations of v^' and w^' with the kernel 

R^ix^xM-A,^) = e qzi - A - m)citl+k( ^ )x , i = 1,2. (77) 

Formulae (1761) transform functions v^l and w^l of three continues variables x, z% and t\ into 
functions (g) (i = 1, 2) depending on continues variable g and discrete variable m. Discrete 
variable m is assotiated with in the kernel of the integral transformation, where m = 
1,...,2N and N is an arbitrary positive integer. Thus, we may state that the solution 
space to the eqs. (l36H3§|) has freedom of two arbitrary functions, ip^'^iq) and i^^'^iq), of one 
continues and one discrete variable, i.e. one has two-dimensional solution space. Using these 
functions we may approximate (at least formally) two initial conditions, for instance, 

(^| tl =o,^| tl =o) - {^%U^\q)), (78) 

and, as a consequence, we are able to solve the initial value problem (I VP) for the system (}36T 
139]) (which describes evolution of and w^ 1 ') with "approximate" initial conditions. We say 
"approximate initial condition" since functions ijr^a (g) (i — 1,2) have one continues and one 
discrete variable, while both variables must be continues in order to represent arbitrary initial 
conditions precisely. We conclude that one provides the same variety of solutions to eqs. ()36H39"j) 
as Sato approach does to the classical S'-integrable PDEs [39] . 

5.2 Solutions to the system (I62U63I) and to its scalar version (CD) 

Eqs. (l47j) have been derived as an algebraic system for the fields v^> and w^> in this case. We 
split the system ( HTj) into two subsystems. First one corresponds to n — 1: 

iVo no 
i=l j=l 

1 = 1,.. .,2N Q , m = 1, . . . ,n . 
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This is the system of 2NqUq linear algebraic M x M matrix equations for the same number 
of matrix fields v^ v,l > and w^' l \ i = 1, . . . , No, j = 1, . . . , no- Namely eqs. (1791) yield functions 
u = and q = w^ 1 ' 1 ' + v^ 1;1+ni ^ as solution to the system (I62|63p . Second subsystem 

corresponds to n > 1 in (1471) : 

No 

qf(l;nm) J^l;m) = [\gfa(n-ni(i))m) _|_ ^,(l;(n-n 2 {i))m) j ^(l;m)^-i+l ^ ^gg^ 

i=0 

= o, if, i < o, 

1 = 1,.. .,2N , n,m = 1, . . . ,n , 

This equation expresses recursively the functions \|>( /;nm ) 5 n > 1, in terms of the functions \|>( Z;lm ) 
and their x-derivatives. The simplest case corresponds to rij(i) = 1, Vz, j (Frobenius reduction). 

Functions v|/( /;nm ) and A^ ;m ^ are solutions to the system (I48H50P with n = 1, 2 (compare with 
eqs. ( 16611671) ): 

^% m \x) = E^r ) (^-^ ;m)t i'^-A^ ;m) t 2 )e fc "r ) ^ (81) 

Af m) = E^ m \z 1 -Af m k 1 ,z 2 -Af m h 2 ), (82) 
1 
2 



c ;2) =^ ^-v(^) 2 + 4 « A r ) + 4 ^ 



m = 1, . . . , n , Z = 1, . . . , 2iVo, a, (3 = 1, . . . , M. 

Here ip^T (yi,y 2 ) and E^ m \yi,y 2 ) are arbitrary scalar functions of two scalar arguments. 

Since the functions \I/(' ;lm ) depend explicitely on the functions A^' m \ the functions v^> 
and w( v 'fi depend explicitely on A^'" 1 ) as well. Functions A^'" 1 ) describe the break of the wave 
profiles because they are solutions of eqs. (l8"2]) . Consequently, the functions v^> and w^> (in 
particular, solutions to the nonlinear PDEs (I62|63p ) exhibit the break of the wave profiles as 
well, unless A^ ;m ^ = const VZ, m. One gets explicite solutions in the later case. 

Consider the scalar case corresponding to eq.flT]), i.e. M = 1 and diagonal matrices a, v, fx 
and A^ ;m ) become scalars. Eqs. fl81l82j) read 

2 

m^ lm )(x) = V> (im %i - A (l ' m) t 1 ,z 2 - A(^)t 2 ) e fc(im;i)x (83) 

i=l 

^(l;m) = E (lm)( Zl _ A&™) tl , z 2 - A ( ' ;m) t 2 ), (84) 

A; ( ' m;1) = l - [v + \J D 2 + Ahk^ m ) + Aj^j , k {lm ' l) = ~ (v - \jv 2 + A~ak^ m ) + 4/2 
m = l,...,n , Z = l,...,2iVo. 
Reduction (j3J) corresponding to the eq.fll]) reduces the eq. (l83l) into the following one: 



2 



i=l 



m = 1, . . . ,n , 1 = 1,.. .,2N , 
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where ip^ 171 ^ and A^'" 1 ) are arbitrary scalar constant parameters. Emphasise, that A^'" 1 ) does 
not depend on variables x and there is no arbitrary functions in the available solution manifold. 

5.2.1 Dimensionality of the solution space 

Let us define restrictions, generated by our algorithm, on the solution space to the nonlinear 
PDEs f l62|63l) . First of all we represent the formal solution to the system flHUl) as follows: 

$ (i;nm) = $( iilm )?< ta ', P {lml) =I M , 1 = 1,..., 2N , n,m = l,...,n , (86) 

reflecting the fact that eqs. ([8"Uj) must express 

we may represent eq. (!79l) in the following compact form: 

JVo n N n 

qj(l;lm) ^(l;m) _ v (i;j)jj(ij;lm) _|_ w (i;j)jj(ij;lm) _ y(r,j)f[(ij;lrn) ^ ^gy^ 

i=l j=l i=l j=l 

1 = 1,.. .,2N , m = 1, . . . ,n , 

where 

i = l,...,N , j,m = l,...,n , I = 1, . . . ,2N . 
Introduce operator f[( ml >3' 1 ) by the following formula: 

2N n 

j2J2fi mm) n {ml ' M) = SiJjkhi, fi (ml ' kn ^ = ( nS mi;ta) fj( m ^™) ) , (89) 

1=1 m=l 

i,n = l,...,N , j,k = l,...,n . 
Then eq. fl8T|) yields: 

2N no 2No n 

v (n-,j) = ^2 ^ i ' im )x( i ' m) fi { ™ hjn \ w {n ' j) = X] # ( ' ;im) A (i;m) n^ mZ;i ™ ) , (90) 

m=l i=l m=l 1=1 

n = l,...,N , j = 1, . . . ,n , 

Now let us turn to eq. flHTj) . We represent arbitrary functions (yi, 1/2) in the following 

integral form 

(1/1,1/2) = | dgi^ 2 ^7 ) (^'^)e 9 ^ 1+92?/2 , (91) 
/ = 1, . . . ,2A^ , m=l, ...,n 0j a, /3 = 1, . . . , M, 

where V'-cT (<Zi 3 ^2) are arbitrary functions of two arguments and one integrates over the whole 
two dimensional space of variables q\ and q 2 which are complex in general. Then eq. flSTl) may 
be written as follows: 

*S lm) ( f ) = [ dq 1 dq 2 ^'' i \q 1 ,q 2 )e giZl+q2Z2 - A ^ (92) 

i=i J 

1 = 1,..., 2N , m = 1, . . . , no, a, /3 = 1, . . . , M. 
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Substituting expression (J92D for ¥aiy lm) into eqs. (190]) we write them in the following form: 

v^ j \x) = (93) 

M 2N n 2 „ 

E E E E / ^g 2 ^ 7 m %i, 52^+^^ 

7 =1 m =l Z=l j=l ^ 

K(3 ( X ) = 
M 2N n 2 „ 

E E E E / d qi d q2 ^<\ qi , q2 )e^^^ 
7=1 m=l 1=1 1=1 J 

n = 1, . . . , JV , j = 1, . . . ,n . 

Now we show that the solution space to the five- dimensional matrix eqs. (I62II63P has two 
arbitrary M x M matrix functions of two continues and one discrete variable (compare with 
Sec J5.1.1 l ) . For this purpose we consider formulae (1931) as integral representations of and 
w ( n <j) with the kernel 

R®(x,z 1 ,z 2 ,ti,t 2 -,q 1 , q2 ,A^) = e ^-i+^ 2 -A^ m) (^i+^ 2 )+A«, ; i = 1,2. (94) 

Formulae (j93l) transform functions u^a and iw^g' 7 '' of five continues variables x, and tj 

(i = 1,2) into functions (51,52) = 1,2) depending on two continues variables 51, 

52 and two discrete variables /,m. However, two discrete variables I and m are assotiated 
with single discrete variable M l]m ^ in the kernel of the integral transformation, where I = 
1, . . . , 2N , m — 1, . . . ,n and N , n are arbitrary positive integers. Thus, we may state that 
the solution space to the eqs. (16211631) (describing the evolution of v^ 1;1 ') has freedom of two 
arbitrary functions, (51, 52) and ip^' 2 \qi, 52), of two continues and one discrete variable, 

i.e. one has three-dimensional solution space with two arbitrary functions of three variable. 
Note that this freedom is not enough in order to solve, for instance, the IVP for the matrix 
system f)62|63p even with "approximate initial conditions" like it was done in Sec l5.1.11 In fact, 
in order to approximate arbitrary initial condition one needs an (infinitely) big number of such 
arbitrary functions of three variables. Thus, further modifications of our algorithm with the 
purpose to increase the dimensionality of the solution space to the derived nonlinear PDEs 
becomes an important problem to be resolved. 

Solution space to the scalar equation ([1]) has the same dimensionality so that formulae 
(I86H94P remains valid with M = 1. 

Consider reduction ([3]) and estimate the dimensionality of the available solution space to 
the (2+l)-dimensional PDE (jlj). As we have seen in both this subsection and Sec l5.1.11 the 
dimensionality of the solution space is completely defined by the kernel of the integral trans- 
formation, see eqs. (j77]l and fl94"|) . Reduction (j3J) reduces eq. (l94"|) with M = 1 to the following 
one 

R {l \x,y,t; Att™>) = ^^^-A^MA^))^ { = 12 (95) 

This kernel has one discrete parameter A^' ;m ^ and, as a consequence, the available solution space 
is one-dimensional with two arbitrary functions of single discrete variable. Number of arbitrary 
functions is predicted by two different values of the superscript i in eq. (l9"oT) . 
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5.2.2 Simplest example of the solution to eq.fll]) 

The simplest nontrivial example of the solution to eq.([T]) is assotiated with N = n = 1. Then 
eg. (185]) reads 

= ^^O^-A^tL^-A^^e* 01 ' *, 1 = 1,2, (96) 
i=i 

where A( /;1 ) (Z = 1,2) are solutions to the nonlinear algebraic equations (18"1"|) which reduce to 
the following ones: 

A (J;i) = E^(z 1 -A^ 1 k 1 ,z 2 -A {l ' 1 h 2 ), 1 = 1,2. (97) 



Then eqs. (1791 become the system of two following equations: 

tf(J;ii) A (J;i) = w (l ; l)^ ; ll) + j = 12 

It has the following solution: 



u = v^ = ^, A 
A 



^(l;ll) 
$ (2;11) ^(2;11) 



(98) 



(99) 



^((-ir +i ^i-(-iy +i ^ 2 ). 



^ = ^V /z>2 + 4 « A( ' ;1) + 4 ^ yi = ^i-A ( ' ;1) *i, i/ 2 = ^-A {i;1) t 2 , 



Ax 



^(l;ll) A (l;l) ^(1;11) 
^(2;11) A (2;1) ^(2;11) 



(A (1;1) _ A (2;l)) e 



•it' 



(i;i) 



A, 



A, 



^(l;ll) A (l;l) 
^(2;11) ^(2;ll) A (2;l) 



(100) 



Formulae (I99|I100I) have six arbitrary functions of two variables ^ (l/i; 2/2) ? A"' 1 ) = E^ ll \yi,y2), 
l,i = 1,2. It is evident that solution u has no singularities if A 7^ 0, i.e., for instance, if 



K 1 >K 2 >0, V (U;2) <0, ^ 11;1 ),V (21;1) ,^ (21;2) >0, Vt/i, 



y 2 - 



(101) 



Since A^ ;1 ^ is implicitly given by the eq. (1971) . constructed function u describes the break of the 
wave profile unless A^ ;1 ) = const. 

Write the explicite form of the solution corresponding to 



Ko = 



A (2;D 



v l + 4/i 
4a 



const, 



(102) 



^ (11;1) (2/i,2/ 2 ) = £1(2/1,2/2) > 0, V> (11;2) (2/1,2/2) = -6(2/1,1/2) < 0. 
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One has 



= (A (1;1) ~ A (2;1) ) (e^'gifoi.lfa) - e- K ^ 2 (yi,V2)) , m\ 

Let us consider the particular case A 1 -'' 1 - 1 = const, 1 — 1,2. Then the solution to eq.(DQ) is 
given explicitely by the eq. (1991) with four arbitrary functions of two arguments ^ ll ' ,l \yi, y 2 ), 
I, i — 1, 2. Consider solution given by formula (11031) with Im(K\) = and Jm(£i) = 0, i — 1, 2. 
Function M| 2i = C onsf,ti=consi is kink. Behaviour of w| x=CO nst as a function of z% and tj depends on 
the shapes of the arbitrary functions £1(3/1,2/2), i = 1,2. In particular, u x=const may be soliton. 

5.2.3 Simplest solutions to eq.((4]) 



In accordance with eq. ([85]) . there is no arbitrary functions in the constructed solution space. 
The simplest nontrivial example of the solution to the eq.([T]) is assotiated with N = n = 1. 
As we have seen above, reduction ([3]) requires 

= const ^ i = i t 2. Eq. (JHSD reads 

= ^(tt;i)/ 1, "M(««y-(A(WFt) | = l!2) ^(ii,-o = oonflt> (104) 
i=i 

Eq. (1981) remains correct. It yields 

„ = ^ (1;1) = ^, (105) 
2 

A _ ^(-l) i+1 iri - (_l)i+ 1 X 2 j e ^(9i/ 2 +92/2)+(-l) l + 1 ^i9i+(-iy+ 1 ^292^(ll;0^(2l;i) ) 

ft (y,t)=x-Aft 1 )y-(A^) 2 t, 

2 

Ai = ^(l;l) _ ^(2;l)^ ^ e £; ^i/2+g2/2)+(-l) l+1 X 1 g 1 +(-l^+ 1 A- 2(?2 ^(ll;i)^(21;i)^ 

It is evident that solution u has no singularities in the case (llOip . 

Write the explicite form of the solution corresponding to conditions H102[) with constant 
% = 1,2: 

_ (A^)-A^))(e^6-e-^£ 2 ) 
If the parameter is real, then the function u is kink. 



6 Conclusions 

The remarkable relations among C-, Ch-, Si- and ^-integrable nonlinear PDEs have been 
found in |2J using a version of the dressing method. We represent a modification of this dressing 
algorithm allowing one to involve a simple example of the Ss-integrable scalar equation (i.e. 
eq.([Tj)) in the above list. In addition, a matrix version of this equation is derived. We have found 
that (at least on the restricted manifold of solutions, as it is described in Secj5]) eq. (pQ) originates 
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from the GL(N S2 , C) SDYM (which is ^-integrable equation) after the special differential 
reduction followed by the Frobenius type reduction imposed on the matrix field. In turn, 
GL(N 2 ,C) SDYM may be derived from the (l+l)-dimensional hierarchy of C/i-integrable 
PDEs [H12]. Thus, the important result is that S3-integrable equation (pQ) has been derived from 
the (l+l)-dimensional hierarchy of nonlinear C/i-integrable PDEs by means of the following 
set of tree reductions imposed on the matrix fields: 

FR — >DR(x) — ► FTR. (107) 

The first and the third reductions in eq. (1107j) introduce the new discrete variables into the 
hierarchy of commuting nonlinear PDEs. The second reduction introduces the new independent 
variable x. After the third reduction, one must to couple two equations of the commuting 
hierarchy of chains in order to write the complete system of nonlinear PDEs. 

Solution spaces to both S^-integrable PDEs with differential reduction and S , 3 -integrable 
system are considered. It is shown that our algorithm allows one to construct a big manifold 
of solutions to eq.([T]) and to its matrix generalization (I62II63P involving solutions with break 
of the wave profile and soliton-like solutions. However, solution space to iS^-integrable PDEs 
available by our method is not full so that IVP may not be formally solved. The problem of 
improving this algorithm with the purpose to describe the full solution space is an important 
problem for further study. 

It is clear that the chain (I107P may be prolonged: 

FTRx(Zi) — ► DRi(xi) — ► FTR 2 (/ 2 ) — ► DR 2 (i 2 ) — ► • • • (108) 
— > DR^-i(xir-i) — > FTR K (l K ). 

Each FTRj(/j) introduces the new discrete parameter l i: while each DRj(xj) introduces the 
new independent variable Xi into the hierarchy of chains of nonlinear PDEs. After the last 
reduction, one has to couple two (or may be more) chains from the hierarchy in order to write 
the complete system of nonlinear PDEs. Study of these PDEs and their possible relations with 
the classical integrable models is postponed for the future work. 

Another evident modification of our algorithm is generalization of the differential reduction. 
Namely, instead of the second order linear PDE ( |22|) one can use i^th order linear PDE with 
arbitrary positive integer K: 

K-l K-j-l 

QKyim) = a {ij) di^ m) {A (m) ) j , m = 1, . . . , 2N . (109) 

Assotiated nonlinear PDEs will be considered elsewhere. 

Author thanks Professor P.M.Santini for useful discussion. The work was supported by the 
RFBR grants 07-01-00446 and 06-01-92053 and by the grant NS-4887.2008.2. 
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